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Abstract 

The dielectric susceptibility of most materials follows a fractional power-law frequency 
dependence that is called the "universal" response. We prove that in the time domain this 
dependence gives differential equations with derivatives and integrals of noninteger order. 
We obtain equations that describe "universal" Curie-von Schweidler and Gauss laws for 
such dielectric materials. These laws are presented by fractional differential equations 
such that the electromagnetic fields in the materials demonstrate "universal" fractional 
damping. The suggested fractional equations are common (universal) to a wide class of 
materials, regardless of the type of physical structure, chemical composition or of the 
nature of the polarization. 
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1 Introduction 



A growing number of dielectric relaxation data show that the classical Debye behavior [TJ [21 [3] 
is hardly ever observed experimentally [U El El E]- Instead it has been derived [U El El H] 
that power laws are a common feature of the dielectric response of most materials for wide 
frequency ranges. The fact that different dielectric spectra are described by the power laws is 
confirmed in many measurements [H El E] for a wide class of various substances. The dielectric 
susceptibility of most materials follows, over extended frequency ranges, a fractional power- 
law frequency dependence, which is called the law of "universal" response [H El E]- This law 
is found both in dipolar materials beyond their loss-peak frequency, and in materials where 
the polarization arises from movements of either ionic or electronic hopping charge carriers. 
These power-law responses are most easily displayed in terms of the dielectric susceptibility 
x(u;) = x'(^) — ix"( u ) as a function of frequency uj. It has been found [HIE] that the frequency 
dependence of the dielectric susceptibility follows a common universal pattern for virtually all 
kinds of materials. The behavior 

X 'H~o; n - 1 , X »~o; n - 1 , (0 < n < 1, cu»^ p ), (1) 

and 

X '(0)-xV)~^ W , *"M~w m , (0<m<l, u<£u p ), (2) 

where x'(0) * s the static polarization and u p the loss-peak frequency, is observed over many 
decades of frequency. Expressions ([T|) and ([2]) serve as the definition of the universal response 
behavior. 

Note that a consequence of the power laws is that the ratio of the imaginary to the real 
component of the susceptibility is independent of frequency. The frequency dependence given 
by equation (Tj[|) implies that the real and imaginary components of the complex susceptibility 
x(u) = x'(w) — iy?'(to) obey at high frequencies the relation 

y"(oj) /7rn\ , 

!^ = cat( T ), («><*). (3) 

The experimental behavior of equation (T5]) leads to a similar frequency-independent rule for 
the low frequency polarization decrement: 

_|lJ_ = tan (_), (4) 

This being a unique consequence of Kramers- Kronig relations and does not depend on any 
particular physical process. 



2 



2 Fractional equations for laws of universal response 

For the region uj 3> uj p , the universal fractional power law ([1]) can be presented in the form 

X(u) = X a (0 < a < 1) (5) 

with some positive constant Xa and a = 1 — n. Using 

{ioj) a = \u\ a exp{i a tt sgn(u>)/2}, 
it is easy to derive relation ([3]). The polarization density can be written as 

P(t, r) = T- x (P(w, r)) = e^- 1 (x(w)E(w, r)) = e Xa ^ ((^)- Q E(^, r)) , (6) 

where P(o>, r) is a Fourier transform JF of P(t, r). Equation ([6]) can be represented by integrals 
of noninteger order a — 1 — n. The fractional Liouville integral [101 EH] is defined by 

The Fourier transform JF of this integral is given (see Theorem 7.1 in [in] and Theorem 2.15 in 
by the relation 

(^i?/)H = 7T^W)(w). 

As a result, the fractional power law (J5]) gives 

P(t, r) = e oXa (I«E)(*, r), (0 < a < 1). (7) 

For the region <C co> p , the universal fractional power law (j2J) can be presented as 

= X(0) - X p (0 < P < 1) (8) 

with some positive constants X/3> x(0), and [3 = m. It is not hard to prove that equation 
is satisfied. The law (jSj) can be presented by the fractional Liouville derivative [TUJ E] that is 
denoted by D+. The differential operator Z?+ of order ft is defined by the equation 

The Fourier transforms T of this derivative (see Theorem 7.1 in [10J and Theorem 2.15 in |llj ) 
is given by 

{FD$f)(u>) = {iuf{Tf){u). 
As a result, the fractional power law (jHJ), gives the polarization density 

P(t,r)=e ^" 1 (xME(^,r)) (10) 
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in the form 

P(t, r) = e x(0) E(t, r) - e oX (3 (^E)(t, r), (0 < (3 < 1). (11) 

Equations © and (|TT|) can be considered as the universal response laws [U El [6] for time- 
domain. These equations allows us to derive fractional equations for electric and magnetic 
fields. 

3 Fractional equations of the Curie-von Schweidler law 

Using ([7]) and ffTTl) . the polarization current density 

, . 9P(t,r) , \ 

J P oi(t,r) = —^ (12) 

can be described by the fractional equations 

J P oi(t,r) = e oXa (£>}T a E)(i,r), (0 < a < 1), (13) 

and 

Jpoj(t, r) = e x(0) DlE(t, r) - e oX p P+ +/3 E)(t, r), (0 < (3 < 1). (14) 

For constant electric field E(t, r) equations ([TBI and (fT4]) show that the time dependence of 
the relaxation of the polarization current density (fT2"I) after the sudden removal of a polarizing 
field follows the power laws, which is widely observed in practice [4] and is known as the Curie- 
von Schweidler law P2J HB]- F° r the changeable field E(t,r), equations (1T3"]) and (114"]) can be 
considered as a generalization of the well-known Curie-von Schweidler law. Let us consider 
some examples of this generalization. 

(1) Using (ITS]) and (fT4"j) . we can derive the usual Curie-von Schweidler law. The most 
elementary of the applied field E(£, r) is the step function 

/ s fo, t < a, 

E(t,r)= . (15) 
[E (r), t>a. 

In this case, equations ([TBI and (]T4"]) give 

3pd(t,r) = e oXa E (r)( a Dl- a l)(t) = 

= e X a E (r) (t ~ a) , (t>a, < a < 1), (16) 
r (a) 



and 



J P oi(t,r) = - W E (r)( a A 1+/3 l)W = -e X/3E (r) v r( _ ; ^ , (t > a, < (3 < I) 
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where we use \TU\ the relation 



( a D?l)(t) = %- o) "° (t>a, a>0). 



Til -a) 



Here a Df is the fractional derivative 



dt k ' 



d 



k r* 



u(t')dt' 



T{k-a)dt k J a (t-t') a - k+v 



where n = [Re(a)] + 1, i.e., n — 1 < a < n. 

As a result, we obtain the usual Curie- von Schweidler law that is described by 

E(t,r) = £ XaE (r) v '- — , (t > a, < a < 1), 

r(a) 

(t- a )- /3 - 1 

E{t,r) = -eoXpMr V , (* > a, < /3 < 1). 

(2) The experimental applied field E(t, r) can be presented as 



E(t,r) = E (r) sin(At). 



Using the relation [10J: 



D° sin(At + 0) = X a sin(At + + twr/2), (a > 0), 
equations (jI3ft. (ITi|l and (fT9]l give 

•W*, r) = e oXa E (r)A a sin(Ai + (1 - a)7r/2), (0 < a < 1), 

and 

•Wt, r) = e E (r) [x(0)A cos(At) - X/3 A 1+/3 sin(At + (1 + /T)vr/2)] , (0 < f3 < 1) 
Equation f[2"Tj) can be rewritten in the form 

J P oi(t, r) = e a(p) E (r) sin (At + 6(/3) 
where a{(3) and describe the amplitude and phase changes by the equations 

■A{py 



Here 



a(P) = ^A 2 (f3) + B 2 (f3), b((3) = arctan 



A(/3) = x(0)A- X/3 A 1+/3 cos(/37r/2), 



B(P) 



B{0) = xA l+P sin(/37r/2). 

(3) For the applied field 

1 0, t < a, , s 

V(t,r)={ (23) 
[E (r)g(t), t>a 

with some function g(t), exact expressions for the polarization current density J po i(t,r) can be 
derived by using the list of fractional derivatives of the function g(t) (see Tables 9.1-9.3 in |10j). 
For g(t) — (t — a) s , where s > — 1, 



r> + i) 

r(s + l-a) 



a D«g(t) = a D?(t - a) s = ^ _\ - '— (t - a) s ~ a , (s > -1, a > 0). (24) 



The fractional derivative of g(t) = cos[A(t — a)] is 

a D?cos\{t-a) = ^ [ x Fi(l, 1 - a, i\(t - a)) + iF x (l, 1 - a, -i\(t - a)) 
21 (1 — a) L 

where i-Fi(a, b, c) is a hypergeometric function [Hj. For g(t) = exp(-Xt), we use 

aJ D a e~ A * = e" A '(t-a)- Q E ljl _ a [-A(t-a)], (0 < a < 1), (25) 
where E a ^[z] is the Mittag-Lemer function [15] : 



fc=0 



If a = (3 — 1, then £1,1 [z] — ex P(^) ; where T(/c + 1) = k\ for positive integer k. 

As a result, fractional relations (fTBl and (Tllft can be considered as a generalization of the 
formulation of the Curie- von Schweidler law from a constant electric field into changeable fields 
E(t,r). 

4 Fractional Gauss's laws for electric field 

Time-domain laws are presented by the fractional integral and differential equations (JTj) and 
( ITTT) . Using the equation 

D(t,r) =e E(t,r) + P(t,r), (27) 

and Gauss's law 

(iiv D(t, r) = p(i, r) 
for the electric displacement field D(t, r), we get 

e dit) E(t, r) + div P(t, r) = p(t, r). (28) 
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Substitution of (ED and (JTTJ) into (ESI) gives 

e Z(t, r) + e X a (i?^)(*, = P(*, r), (0 < a < 1), (29) 

e„x(0)£(t, r) - e QX( 3(D p + Z)(t, r) = p(t, r), (0 < < 1), (30) 
where Z(t,r) = divE(t,r). Using (D+ I+Z)(t,r) = Z(t,r) (see Lemma 2.20 in [11]), we obtain 

(D«Z)(t,r) + X *Z(t,r) = -(D«p)(t, r), (0 < a < 1), (31) 

(pJZ)(f,r)-^Z(f,r) = -— p(f,r), (0 < /3 < 1), (32) 

which are fractional differential equations of Gauss's law for the electric field E(t,r). 
For a fixed (stationary) region R of medium, we define the total electric charge 



Q(t) = / p(*,r)<fV. (33) 

J R 

The electric field E = E(t, r) passing through a surface S = <9i? gives the electric flux 

$ E (t)= /(E,dS)= / dwEdV. (34) 



The integration of equations fl3Tl) and fl32|) over the region i? gives the fractional equations 

(£>?$s)(t) + xA(t) = (0 < a < 1), (35) 

(L>?<^)(t) - ^S*(t) = — (0 < P < !)• (36) 

These equations represent the integral Gauss's laws for the electric field in dielectric media. 
Note that .D? is the differential operator of order a that is defined by equation Q. If E(£, r) 
is defined by 

fo, t < a, , s 

E t,r = < ' ' 37 

\E(t,r), f>a, 

then the operator transforms into the operator a .D" that is defined by equation (fT7|) . 

Consider equations fl35l) and fl36|) in the form 

£>Xt) - Ait(t) = /(*), < a < 1, (38) 

where presents Z(t,r) or A is — Xa or x(0)/X/3> the function f(t) is (1/ffo) a.Dfp{t,r) 

and (— 1/eoXp) p(t, r )> or (l/^o) aDfQif) and (— l/eaXp)Q{p)- The fractional derivative a .Df is 
defined by (JTZD- Note that 

(£>»(*) = lim ( a D?u)(t). 

a— >— oo 
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J.H. Barrett [16] in 1954 first considered the Cauchy type problem for a linear differential 
equation f[38l with the initial conditions 



( a Dr L u)(a)=C } (39) 

where a D^~ l = J-\~ a is the Riemann-Liouville fractional integral [llj. Note that any change 
in the past of the input function of a fractional order system changes the future of the solution 
[T7J [HJ US]- Hence the past of such systems cannot be represented by a finite set of local 
conditions [T7J HH [19] , and the initial conditions have the integral form (1391) . 

Barrett proved in [16J (see also [11] Theorem 4.1 and Example 4.1.) that if f(t) is an 
integrable function on (a, b), then the problem (I3"5j) . (13"91) has the unique solution given by 

u (t) = C(t - a) a - l E^ a [X{t - a) a ] + f (t - t') a ~ l E a , a [\(t - t') a ] /(*') dt', (40) 

J a 

where E a ^[z\ is the Mittag-Leffler function defined by (|2"6"1) . For f(t) = 0, we obtain 

u(t) = C(t - aY^E^l^t - a) a }. (41) 



To consider the asymptotic behavior of the solutions (140]) and fl4T|) . we can use the integral 
representation [201 ED E2] of the Mittag-Leffler function: 



i r e ?1/< Y (1 -' 9)/Q 

E ^ = 0^ / / y «*6, (1 < « < 2 )' (42) 

where 7ra:/2 < 5 < min{iT,7ra}. The contour 7(0, 5) consists of two rays S 1 ^ = {arg(£) = 
—5, |£| > a} and S+,5 = {arg(£) = +5, |£| > a}, and a circular arc C$ = {|£| = 1, — 5 < arc(£) < 
5}. Let us denote the region on the left from 7(0, 6) by G~(a, 5). Then the asymptote of (l4"2l 
has the form [2T1I22] : 

OO _t 

= -X) r(/-aA;) a (kl ^ 00), (43) 

and 5 < | arg(z)| < 7r. In our case, z = X(t — a) a , arg(z) = tt. As a result, we arrive at the 
asymptote of the solution, which exhibits power-like tails. These power-like tails are the most 
important effect of the noninteger derivative in the fractional equations. 



5 Conclusion 

In this paper, it has been shown that the electromagnetic fields in a wide class of dielectric 
materials must be described by differential equations of fractional order with respect to time. 
The orders of these equations are defined by exponents of the "universal" response laws for 
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frequency dependence of the dielectric susceptibility. A remarkable property of the dynamics 
described by the fractional equations for electromagnetic fields is that the solutions have power- 
like tails. The typical features of the "universal" electromagnetic phenomenon and the suggested 
fractional equations are common to a wide class of materials, regardless of the type of physical 
structure, chemical composition or of the nature of the polarizing species, whether dipoles, 
electrons or ions. 

For small fractionality a (or /3), it is possible to use a e-expansion [23j over the small 
parameter. The suggested fractional differential equations, which describe the fields in dielectric 
media with power-law frequency-domain response, can be solved numerically. There are several 
numerical methods to solve fractional equations (for example, see PHI [2U [221 ES] ) • In this paper 
we consider only general equations without numerical examples. 

The presented approach does not take into account the dielectric dispersion region around 
the loss-peak frequency uj p . This is a limitation of the fractional calculus formalism. Note 
that the universal dispersion is only valid over a finite frequency region. There exist numerical 
methods without a piori assumptions to tackle frequency to time and time to frequency con- 
versions using distribution of relaxation times approach. There are several papers on this topic 
by different authors (for example, see (23 123 123 EDI EH)- 

Note that Tuncer has shown [32] that in binary mixtures Maxwell- Wagner polarization could 
lead to a universal dispersion on a finite frequency region. In that study the disorder in the 
system leads to a universal dispersion through the Gauss' law. 

Note that it will be interesting to find a generalized physical explanation for fractional 
power laws. Unfortunately, after more than 30 years we still do not have an explanation of 
the universal response from first principles. There are several papers on the fractional kinetic 
approach to the dielectric universal response by different authors (for example, see [3"3"|I3"4"1 I3"5]). 
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